This paper, based on panel data of 16 cities of 2007, 2012 and 2017 in Shandong province, establishes a coordination model of population-economy-environment based on entropy and physical coupling model, which is used to study the temporal-spatial evolution of population-economy-environment system coupling coordination and to estimate their spatial Moran's indexes, and presents suggestions for Shandong's future development. Among 16 cities in Shandong province, their coupling coordination is improving thanks to a rising comprehensive development. The western inland is far behind the eastern and central in economy and population. A low spatial autocorrelation exists in their economy-population-environment development, but weakly positive between comprehensive development and coordination. Regional coordinated potential needs to be further developed.